Factorizing the time evolution operator 

P. C. Garcia Quijas and L. M. Arevalo Aguilar* 
Centro de Investigaciones en Optica, Loma del Bosque No. 115, 
Fracc. Lomas del Campestre, Leon, Gto. Mexico. 

Abstract 

\& ' There is a widespread belief in the quantum physical community, and in textbooks used to teach 

o 

Quantum Mechanics, that it is a difficult task to apply the time evolution operator e itH / h on an 
initial wave function. That is to say, because the hamiltonian operator generally is the sum of 
two operators, then it is a difficult task to apply the time evolution operator on an initial wave 
function t/)(x, 0), for it implies to apply terms like (a + b) n . A possible solution of this problem 
(f-) ■ is to factorize the time evolution operator and then apply successively the individual exponential 

>: 

operator on the initial wave function. However, the exponential operator does not directly factorize, 
i. e. e a+b ^ e a e b . In this work we present a useful procedure for factorizing the time evolution 
operator when the argument of the exponential is a sum of two operators, which obey specific 
commutation relations. Then, we apply the exponential operator as an evolution operator for the 
Qh! case of elementary unidimensional potentials, like the particle subject to a constant force and the 

harmonic oscillator. Also, we argue about an apparent paradox concerning the time evolution 

d \ 

^ { operator and non-spreading wave packets addressed previously in the literature. 
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I. INTRODUCTION 



Quantum Mechanics is a successful theory. Although highly counterintuitive, using it 
we are able to explain the microscopic world. Also, Quantum Mechanics have discovered 
many natural process which have culminated in practical technological applications, like the 
transistor and the Quantum Cryptography. However, Quantum Mechanics is a difficult field 
of study. For it takes many years to develop the necessary skills to understand their relevant 
concepts. 

One of the hardest skills to develop is to understand the technique to solve the fundamen- 
tal equation of Quantum Mechanics, i. e. the Schrodinger equation. In fact, there are few 
cases where this equation was analytically solved. One of the principal factors that impede 
the straightforward solution is that this equation involves not usual mathematical concepts, 
like operators: 

iftjU(t) = H*(t), (1) 

where the Hamiltonian H is an operator that has to be self-adjoint, and it is generally the 
sum of two or more operators, let us say: 

H = a + b. (2) 

In the literature, the most used technique to solve Equation (1) is to find the eigenvalues 
and eigenf unctions of the time independent Schrodinger equation: 

Htp n (x) = E n ip n (x), (3) 

where E n and if) n (x) are, respectively, the eigenvalues and eigenfunctions of the hamilto- 
nian H [1, 2, 3, 4]. Then, the time dependent wave function is constructed taking the 
superposition of the eigenfunctions of the hamiltonian: 

oo 

V(x,t) = J2^ lEnt/h CnMx), (4) 

n=0 

where c n = f ^/(x,0)ip n dx is the scalar product between the initial state of the system and 
the eigenfuctions of the hamiltonian. In this paper, we call this method the eigenstates 
method. 
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A second way for solving Eq. (1), if we consider a time independent Hamiltonian, is to 
integrate Eq. (1) with respect to time, to obtain [1]: 

t) = e-%H(x, 0) = e A+ H(x, 0), (5) 

where ^(x, 0) is the initial wave vector, A = — (it/h)a and B = —{it/K)b. In this paper, we 
call this method the evolution operator method. 

Essentially, both ways for solving the Schrodinger equation are the same. This can 
be proved by expanding *f?(x, 0) in terms of the eigenfunctions of the Hamiltonian, i.e. 
ty(x, 0) = c n ip n (x), and inserting it on the right hand side of Eq. (5) to produce Equation 
(4). An alternative way to solve the Schrodinger equation is the technique developed by 
Feynman, called the Feynman propagator method [5, 6, 7]. 

The trouble with Equation (5) is that, in general, A and B do not commute. This makes 
difficult to apply the time evolution operator to the initial state vector given in Eq. (5). In 
fact, the problem is how to make the expansion of a function of noncommuting operators 
like that in Eq. (5), i. e. e A+ ^ = E^LoC 1 / 71 ')^ + B) n , in such a way that all the B 
precede the A, or viceversa. This problem has been already studied by many authors, and 
some theorems have been proved to handle this expansion. For example, Kumar proved the 
following expansion for a function of noncommuting operators [8]: 

f(A + B) = jr^C n (A,B)f^(A), (6) 

77/. 

n=0 

where C n (A, B) is a coefficient operator given in terms of B and the commutator [A, B] [8]. 

Also, Cohen has proved the following expansion theorem for the operators x and p [9]: 
Given a function F(x, p) then 

x POO 

F n (x,p) = J2a n k u k (x) / ui(x + 6)e^d6, (7) 

n=0 J -°° 

where a k and u k (x) are the eigenvalue and eigenfunction of the eigenvalue problem 
F(x,p)Uh(x) = akUk(x). In particular, the expansion for the function (Xx + p) n has given 
as [9]: 

(« +»*-££ jkZm W2)*v~*~«-»-y. (8) 

fc=0 s=0 >' \ s I 
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In general, these expansion theorems have produced a high cumbersome expressions that 
are very difficult to apply. 

One of the possible paths to avoid the expansion of functions of two noncommuting 
operators, in the case of the exponential operators, is to factorize the argument of the 
exponential. This approach facilitates the application of the exponential operator because 
now, when the exponential operator is factorized, we have only to expand the exponential of 
a single operator, i. e. e A e B , which is more simple. However, the factorization of exponential 
operators is not an easy task. To our best knowledge, the evolution operator method has 
been applied to unidimensional problems in only four other related articles [10, 11, 12, 13]. 

The main goal of this paper is twofold, first we shall show a procedure to factorize the 
exponential operator and, secondly, we shall show how to apply the factorized exponential 
operator on an initial wave function. The method of factorization that we shall present in 
this paper has been used in Quantum Optics. Also, this method has been proposed as a 
possible tool to improve some misconceptions in the teaching of Quantum Mechanics [14]. 
Therefore, an important objective of this paper is to review this method in order that it 
becomes available for the people outside these fields. 

Although the three methods for solving the Schrodinger equation mentioned above have 
to give the same result, the evolution operator method is, in some way, quite different of the 
eigenstates method and the Feynman propagator. For example, in the eigenstates method 
we need to look for the eigenfunctions of the hamiltonian where the particle is placed; 
on the contrary, the evolution operator method does not give any information about the 
eigenfunctions of the hamiltonian. Also, the Feynman propagator method need to look for 
all the possible paths the particle can take from an initial wave function to a final one, and 
the evolution operator method does not inquire for these possible paths. In some sense, the 
evolution operator method is more direct than the other two. 

In summary, this paper address the problem of how to factorize the exponential of a sum 
of operators, in order to be able to apply it as an evolution operator, when the operators 
obey certain commutation rules. To make the factorization we use the tool of the differential 
equation method [15, 16], which requires that both sides of an equation satisfy the same 
first-order differential equation and the same initial condition. For a review of these tools 
see the work of Wilcox [15] and Lutzky [17]. This method has been used successfully in the 
field of Quantum Optics [18, 19, 20]. We shall show that this method is useful and easy to 



4 



apply in the unidimensional problems of Quantum Mechanics. 

This paper is organized as follows: In Section II we will present the method and show how 
to apply it for factorizing an exponential operator. In Section III we give a specific example 
when the operators obey certain commutation rules. In Subsection A of this section, we 
apply the found factorization to the case when the particle is subjected to a constant force. 
In Section IV we present the factorization of the exponential operator when its argument 
obeys a more complex commutation rules; in subsection A of this section the factorization 
found is applied to the harmonic oscillator. In subsection B of section IV, we derive another 
way to factorize the harmonic oscillator and we show that both factorizations give the same 
evolution function (In Appendix A we derive yet another way to factorize the harmonic 
oscillator). In Section V we address a supposed limitation of the evolution operator method, 
we demonstrate that the limitation is because the initial wave function used to show the 
apparent paradox is outside of the domain of the hamiltonian operator. 

II. THE METHOD 

As the global purpose of this paper is pedagogical, in this section we show how the method 
works. Our intention is that this method can be used for researchers of any field to find the 
evolution state from an initial wave function. In order to be explicit we separate the method 
in three steps and apply it to obtain the well know Baker-Cambell-Hausdorff formula: 



where 5 is a constant. Notice that after we have solved this easy problem we will progressively 
increase the difficulty of the commutation relation. 

To make the factorization of the exponential of the sum of two operators we proceed as 
follows: 

1. Firstly, we define an auxiliary function in terms of the exponential of the sum of two 
operators, its commutator and an auxiliary parameter £: 



e (A + B) 



= e A e B e f ^ 



(9) 



When [A, [A, B}} = [B, [A, B]] = i.e. [A, B] = 5. 



F(0 = e^ A+B \ 



(10) 
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F(0 = e h(0[A ' B] e MoA e m)B . (11) 

Note that in Equations (10) and (11) we are defining separately F(£) as a function 
and its factorization. That is to say, if e fl ^ A,B ^e^ 2 ^ A e^^ B is the factorization of 
e^ A+B \ then F(f) = e ^ A+B) = e fl ^ A ' B k f2 ^ A e fs ^ B . 

2. Secondly, we differentiate Equations (10) and (11) with respect to the parameter £ to 
obtain: 

^f- = (A + B)e^ A + B ) = (A + B)F(£), (12) 



dF (^ = df 1 ^) \\, B]c fl( ^ A ' B] c- m)A c m)B + 

After that, we need to put in order the operators of Equation (13). In order to make 
the new arrange we use the fact that the operators are self-adjoints, i. e. e^ A B = 
(e^Be-^ e^ A , and the well know relation: e^ A Be~^ A = B + i[A, B] + £[A, [A, B]] + 
. . . , see reference [16] . That is, we have to pass the exponentials to the right in the 
right hand side of Equation (13). In this case: 



e M()MAe-^M = A (14) 

and 

e f2(()Afi e -Mt)A = B + 5f ^^ (15) 

where we have used the fact that [A, B\ = 5. If we substitute these relations in Eq. 
(13) we obtain: 

dF ^ = d f^) \\, B]c fl(e)M c h(e)A c f3(e)B + 
di dS, 



d h(0 A c .h(Z)\AB] 



e .f2(0A e M0B + 



d£ 

d Ml e h(OlA,B][ B + 5f 2 (0]e f ^ )A e^ B . (16) 
«4 
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Now, using the relation e h ^ A ^Be~ h{ -^ A ^ = B in Eq. (16) we obtain: 

d f 2 ^ ^ c /i (fl M c h (Z)A c f 3 (£)B + 
Ml [fl ^/ 2 (()] e A©^] e M)V'^. (17) 

That is, we successfully passed all the exponential to the right and we can write Eq. 
(17) as: 

^(0_f4fi(0 A ^ + ^ + ^ + , /a(0] 



Xe Mt)[AA e MQA e MQB m (18) 
That is, by Equation (11) we arrive to the following result: 

dF(o _ !m\ lAB]+ *m A+ 'm B+Sfm \ m . (19) 



3. Finally, as a third step, we must to compare the coefficients of Eq. (12) and Eq. (19), 
from which a set of differential equations is obtained: 



6 ^r +6 ^r m) = °' (20) 

subjected to the initial condition /i(0) = /2(C)) = /3(C)) = 0. In this case the solutions 
are: 

/ 2 (0 = 
/ 3 (0 = e, 

/i(0 = -^ 2 - (21) 
After substituting Eq. (21) in Eq. (11) we arrive to the following Equation: 

e «A+B) = e -^[i,B]/2 e a e ^. (22) 

Setting £ = 1 we obtain the usual Baker-Campbell-Hausdorff formula. 
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This method facilitates the application of the exponential operator, because now we have 
to handle only individual operators function. The proposed factorization of Equation (11) 
is one of the possibilities, also we can define as: 

F(0 = e^V^V 3 ^, (23) 

or make another arrange of the exponentials, as for example e /ll( ^ A e /l2< ^' B e ?i:i< ^ A . Notice 
that Equation (23), in contrast to Equation (11), does not use the commutator in the 
exponential functions. This arrangement could be used to treat specific problems, as the 
harmonic oscillator. 

In fact, when the method is dominated this arrangement is a set of crafted directions, 
which gives a factorization of the evolution operator. In the majority of different 
arrange will produce a different set of differential equations and, obviously, a different set 
of solutions. We give an explicit example of this fact in the case of the harmonic oscillator, 
see Equations (56), (48) and (Al). It is very important not to confuse the Baker-Campbell- 
Hausdorff formula with the method addressed here. Each one represents a different way to 
factorize exponential operators. 

On the other hand, this method can be used to improve some misconceptions of Quan- 
tum Mechanics [14, 21]. Immediately we present the different cases that appear when the 
operators obey different commutation rules. 

III. CASE 1: [A, B] = C, [A, C] = AND [C, B] = -k 

This section is organize as follows: Firstly, we make the factorization of the exponential 
operator when the operators obey the commutation relations given by Equation (24). Sec- 
ondly, in Subsection A we use the factorized exponential to solve the problem of a particle 
subjected to a constant force. 

Therefore, we begin the factorization of exponential operators by analyzing the case when 

[A, B] = C, [A, C] = 0, and [C, B] = -k, (24) 

where A, B and C are operators and A; is a c-number (in general, we use the simbol " to 
denote operators). 
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In the present case, we propose the factorization function as: 

= e^ A+ ^ = e K)A e g ®*e h ® d e r ®, (25) 

by differentiating Eq. (25) with respect to £ we obtain for the left hand side 

dF(0 



and, for the right hand side 



(A + B)F(0, (26) 



dF(0 



di 

where we have applied the fact that 



df j dg ~ dg dh A 



^(0, (27) 



e^Be^ A = B + £[A, B] + [i, S]] + . . . , (28) 

and we have used the commutation relation of Eq. (24). 

By equating the coefficients of Eq. (26) and Eq. (27) we obtain the following system of 
differential equations: 

df(Q . 
di A ' 

subjected to the initial condition F(0) = 0, which implies: 

/(0) = 0(0) = MO) = r(0) = 0. (30) 
By solving Eq. (29) with the initial condition stated in Eq. (30), we finally obtain 

e i{A+B) = e (emk e a e zB e -(emc. (31) 

Setting £ = 1 we obtain the factorization we were looking for. 
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A. Application: A particle subject to a constant force 



One application of the evolution operator method is when we study the time dependence 
of a quantum state. There have been some results in this approach when the operator is 
the energy of a free particle [10], or the energy of a particle subjected to a constant force, 
that is V(x) = —Fx [11]. In this subsection, we use the factorization found above to solve 
the problem of a particle subject to a constant force, whit the help of the Blinder's method 
[10]. The Blinder's method show how to apply the evolution operator like and infinite sum 
for a free particle: 

my i 



g 2mh 



n=0 



^ \ 9,m. I n\ 



d 2 



(32) 



2m J n\ \dx 2 

For a free particle the wave function at time t is obtained by operating with the evolution 
operator on the initial wave function. Taking as an initial wave function: 

1 



0) 



-e w 



(33) 



(^v^) 1/2 

where a is the width of the wave packet. The Blinder's method consist in the application of 
the identity [10]: {z~ 1/2 exp{~x 2 /Az}} = §- z {z~ 1/2 exp{-x 2 /Az}}. This identity allows 
us to apply the evolution operator on initial wave functions like gaussian wave packets, for 
details see reference [10] . 

For a particle subject to a constant force, i.e. V(x) = —Fx, the wave function at time t 
is given by: 



^(x, t) = exp 



it 











*M), 



(34) 



defining A = p 2 (2m) 1 , B — —Fx, and using the commutation relations between p and x 
we can deduce the following commutation rules: 



[A,B] 



ihF 



-p, 



[A, C] = 0, and [C, B] 



h 2 F 2 



(35) 



m m 

where C = ihFp/m. If we identify k = h 2 F 2 /m, then the commutation relations of Eq. 
(35) are similar to that of Eq. (24). Therefore, if we use Eq. (31), we can write Eq. (34) 
as: 



ty(x, t) = exp 



it 3 F 2 
3hm 



exp 



it 



2mh 



P 



exp 



itF 



-X 



exp 



it 2 F 
2mh 



P 



*(x,0). 



(36) 
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Using the theorem exp[£A\F(B) exp[— £A] = F(exp[C,A}B exp[— £A]) [16], we can rear- 
range Eq. (36) as follows: 



ty(x, t) = exp 



iF 2 t 3 
6mh 



exp 



itF 

IF 



-X 



exp 



it 



2mh 



P 



exp 



it 2 F 
2mh 



P 



*(x,0). (37) 



Taking as an initial state that of Eq. (33) we finally obtain: 



/ x -1/2 , 

= crV27r) (l + (iht/2ma 2 )) exp 



(x - (Ft 2 /6m)) 



x exp 



(x- (Ft 2 /2m)Y 



4 (cr 2 + (iht/2m)) 

Equations (37) and (38) are exactly the same equations obtained by Robinett [11]. 



(38) 



IV. CASE 2: [A, B] = C, [A, C] = 2 7 i AND [B, C] = -2 7 5 

In this section, we carry out the factorization of the exponential operator when the 
commutation rules are given by Equation (39). Then, we will show in Subsection A that 
these commutation relations are the same of the harmonic oscillator. On the other hand, in 
Subsection B we show an alternative way of factorization for this problem, and show that 
the evolution given by the evolution operator are the same in both cases. 



[A, B\ = C, [A, C\ = 2 1 A [B, C\ = -2 7 5. (39) 
In this case, using an arrangement similar to Equation (23), we define the function as: 

F(0 = e^ A+ ^ = e /«) V®V<«f (40) 



By differentiating Equation (40) with respect to £, we obtain: 

dF(0 



d£ 



= (A + B)F(t), 



and 



dm 



,(0^-/(0^#- 7^(0/(0 dm 



dg(0 , dh(0 



d£ 



+ 



di 

27/(0^(0 



+ 



di 

dm 

di 

dg(0 



+ 7 2 / 2 (0<? 2 (0 



di 

dm 

di 



c 

B 



di 



+ 7/(0 



dm 

di 



(41) 



A F(0, (42) 
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where we have applied the relation e^ A Be^ A = B + £[A,B] + [A, B]] + . . . . 

By equating Equations (41) and (42), we obtain the following system of differential equa- 
tions: 



^ + ifiO^ + ^ - 27/(?) 9 (0^ + 7 2 /WK)^ = 1, 

subjected to the initial condition F(0) = 1, which means: g(0) = /(0) = /i(0) = 0. By 
solving Equation (43), with the initial conditions, we obtain the following solutions: 



/(0 = M0 = ^tanfo/7/2), (44) 
5 (0 = ^= sin (eVT) . (45) 

Setting £ = 1 we obtain the factorization we were looking for, that is Eq. (40). 

As it was stated at the end of Section 2, the factorization given in Equation (40) is only 
one of many possibilities. Since the operators do not commute, various orderings on the right 
hand side of Equation (40) represent different substituting schemes as we will show in the fol- 
lowing subsections and in the appendix. For example, we can propose a different arrangement 
= e fl ® A e fa ®*e f3 ® A , or inclusive add the commutator C: F(f) = e h ® d e f ® A e 9 ®B. 



A. Application: The one-dimensional harmonic oscillator. 

One of the most important systems in Quantum Mechanics is the harmonic oscillator. For 
it serves both to model many physical systems occurring in nature and to show the analytical 
solution of the Schrodinger equation. The Schrodinger equation for this system has been 
solved in two ways, firstly by analytically solving the eigenvalue equation and, secondly, by 
defining the creation and annihilation operators [1, 2]. We solve here the problem using the 
evolution operator method. This method allows us to find the evolution for the harmonic 
oscillator and avoids to deal with the stationary states. 
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For the one-dimensional harmonic oscillator the wave function at time t is given by: 

it ( p 2 1 



^(x, t) = exp 



H — mw 2 x 2 
h \2m 2 



*(x,0). 



(46) 



Defining A = — (itmw 2 / '2H) x 2 , B = — (it/2mh)p 2 , and using the commutation rules be- 
tween xypwe can deduce the following commutation rules: 



A,B =C, A,C 



-2w 2 t 2 A and 



B,C 



2wH 2 B, 



(47) 



where C 



iw 2 t 2 
2H 



(xp + px). If we identify 7 = w t , then these commutation relations 



correspond to that of Equation (39). Therefore, using the factorization found in Eq. (40), 
the Equation (46) becomes: 



(48) 



where 



ih 

2moj 



tan(a;t/2), 



. imuj . . . 
bit) = — sin(a;t). 



(49) 



For the one-dimensional harmonic oscillator the wave function at time t is obtained by 
operating with the evolution operator, i. e. Eq. (48), on the initial wave function. Taking 
as an initial wave function: 



0) 



- e 4?*, 



(a^y/ 2 

we finally obtain the state of the system at any time t is: 

1 1 



(a^) 1 / 2 y/cos(ujt) + (2ih/mu)(l/4:a 2 ) sm(ut) 



x exp 



{imuj/2h) sin(wt) - (l/2a 2 ) sin 2 (cut/2) + (1/4(T 2 



-x 



cos(wt) + (2ih/muj)(l/Aa 2 ) sin(cjt) 
From Equation (50) we can calculate the probability distribution function: 

1 



(50) 



|tt(x,f)|' 



x exp 



(Jv^F) [cos 2 (cjt) + (2fr/mcj) 2 (l/4CT 2 ) 2 sin 2 (^)] 
2(l/4(7 2 ) 



1/2 



X 



^cos 2 (cut) + (2^/mw) 2 (l/4cr 2 ) 2 sin 2 (cjt) 
In the preceding case we have used the following trigonometric identities: 1 
2sin 2 (wt/2) = cos(cjt) and sin(atf) = 2sin(wt/2) cos(o>t/2). 



(51) 
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B. Another way to factorize the harmonic oscillator 



In this subsection we present another way to factorize the evolution operator the harmonic 
oscillator. Then, we apply the new factorization on an initial wave function. 
In this case we propose the factorization function as: 

F(0 = e^ A+ ^ = e^V^V 3 ^. (52) 

Applying the method of factorization, we obtain the following set of differential equation: 

^ + T/?(0^ + ^-2T/ 1 tt)/ 2 «)^ + 7 2 /f(f)/ 2 2 (f)^f = 1, 

MQ^-fM^-iMOfm^ = o, (53) 

subjected to the initial condition F(0) = 1, which means: /i(0) = /2(0) = fs(0) = 0. This 
set of differential equations is identical to that of Eq. (43). By solving Equation (53), we 
obtain the following solutions: 



/i(0 = / 3 (0 = ^tan(eV7/2), (54) 
h (0 = ^= sin (eVT) • (55) 



Application 

Now, using the factorization given by Eq. (52) to solve the harmonic oscillator problem, 
we obtain the following evolution function: 

= e-^V^e" ^* (a, 0), (56) 

where 

<x(t) = — tanM/2), 

/?(*) = ^sin M . (57) 
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FIG. 1: A plot of the function |^(x,i)| 2 , given by Equation (51). 



If we use again the initial wave function: 



V(x, 0) 



1 




we obtain the following evolving wave function: 



x exp — 




1/2 



(58) 



Equation (58) is exactly the same wave function found in Subsection A, i. e. Equation 
(50). Therefore we can conclude with one of the main points of this paper: The factor- 
ization could be made of different ways and all of them have to give the same 
result when they are applied to an initial wave function. 

The evolution operator of the Equation (56) was also proved by Beauregard [22]. However 
in this work was not used the factorization method and the solution is given by an ansatz. 

In Fig. (1) we have plotted the probability density given by Eq. (51). Were we have set 
er = ^h/muj = 1 and x = 1. 
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V. A NOTE 



Holstein and Swift published a paper in which they presented a cautionary note about the 
usefulness of the evolution operator method for obtaining the wave function at any future 
time t from the one at t — [23]. Notably, in this paper Holstein and Swift showed a par- 
ticular case where the evolution operator method does not works, but if this case is analyzed 
by the eigenstates method it works very well. That is to say, the results obtained with both 
methods do not coincide. Therefore, a contradiction between the evolution operator method 
and the eigenstates method arise, the goal in this section is to present a solution of this 
problem. 

Firstly, we recall the arguments of reference [23]. In their argumentation, they considered 
a "free particle" represented by a one-dimensional wave packet described by the function 

iPh(x,0) = exp [—a 2 / (a 2 — a; 2 )] , for \x\ < a 

= 0, for \x\ > a. (59) 

They argued that iI>h(x, 0) is a "good" function because ?Ph(x, 0) and all its derivatives exist, 
are continuous for all x, and vanish faster than any power as \x\ — > oo. When they apply 
the evolution operator exp ^— j-H^j to the function iI>h(x, 0), they found that 

n=0 v 7 v 7 

when \x\ > a since ipn{x, 0) and all its derivatives vanish for \x\ > a. From this result 
their conclusion was that the particle described by the function iI>h(x, 0) is confined within 
—a < x < a for all time. That is to say, the wave packet does not spread. However, if the 
problem is solved using the eigenstates method then the wave packet do spread, see reference 
[23]. 

In the next two subsections we analyze this argument from two points of view. In sub- 
section A we analyze it from the mathematical point of view. In the subsection B we give 
a physical argument. 

A. Mathematical view 

Mathematically, the argument given by Holstein and Swift is well established. From 
the mathematical point of view, they correctly stated that the function is a well behaved 
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function. That is, because outside of the interval \x\ > a it vanishes, this function does not 
have any singularity and, then, it is an analytical function. Therefore, the conclusion is that 
it is not possible to apply the series of e~ ltH / h to the function tp H (x). In fact, there is a large 
set of such function, i. e. C°°, see reference [24]. 

From this conclusion we can deduce that the evolution operator method fail when one 
apply it on tpn (x). This is a very subtle problem. As it was stated in the previous paragraph, 
the function ipnix) is a well behaved function from the mathematical point of view. There- 
fore, it seems as if the evolution operator method fail when it is applied to an analytical 
function. 

However, after a careful analysis, the only that can be concluded is that the evolution 
operator method fail for non analytical function. For example, the analysis of J. R. Klein 
concludes that the evolution operator method hold for ip(x, 0) in a suitable dense subset, see 
section IV of reference [24]. Another possible argument, stated in the following subsection, 
involves the differences between the Hermitian and self adjoint operators [24] and the fact 
that the particle is free. 

Before to give a possible argument, let us recall a similar function studied by Araujo, et. 
al. [25]. Araujo, et. al., exemplify with the function: 

Uab = e l/(*-«)(*-6) > f or a<x< b, 

u ab = 0, for < x < a and x > b. (61) 

For the entire interval, — oo < x < oo, this function is very similar to the function of Holstein 
and Swift, iPh{x,0). However, Araujo, et. al., use the function (61) only to show that the 
Hamiltonian is not a self-adjoin operator [25]. 

B. Physical view 

In this subsection, we will show that the function used by Holstein and Swift is not a 
valid physical function in the case of the free particle. A crucial point is that unbounded 
operators can not be defined on all functions of the Hilbert space [24, 25]. In order to be 
able to argue this point, we write the following two explicit assumptions given by the cited 
authors [23]: 

1) The particle is a free particle. 
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FIG. 2: A plot of the function iPh{x, 0). 

2) The state of the particle at the initial time is given by the function ipnix, 0). 

Our principal point will be that the statements 1) and 2) can not be true at the same time. 

We begin analyzing the function from the mathematical point of view in the entire domain 
of x, that is {x e R\ — oo < x < oo}, see Fig. (2), where the particle can be found. Notice 
that the function of Holstein and Swift is defined as zero outside the interval \x\ > a, 
therefore, strictly speaking the Fig. (2) is not ipH^x). However, we use the entire interval 
— oo to +oo because the analysis of reference [23] is for a free particle. 

This is a very peculiar function because it has two singularities and many limits: 

Lim x ^ +OQ %lj H (x,{)) = 1, 
Lim x ^ a +^ H (x, 0) = oo, 
Lim x ^ a -ip H (x,0) = 0, 
Lim x ^_ a +tp H (x,0) = 0, 

Lim x __ a -^ H (x,0) = oo. (62) 

From the previous equation we can say that the function is not continuous for all x (remem- 
ber that in this subsection we are studying it in the whole real axe). Furthermore, it does 
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not vanish faster than any power of x as |x| — > oo. Therefore, the function is not continuous 
and \ifj(x)\ 2 dx is not finite, that is, it is not a square integrable function in the interval 
— oo < x < oo. It is important to stress that we are analyzing the case in which the function 
has a definite value outside the interval \x\ > a. 

Now, Holstein and Swift [23] define the function as zero in the interval \x\ > a. With 
this restriction the function is mathematically a "good" function and can pass as a valid 
function. That is, as it was stated in the previous subsection, in this case the function does 
not have the behavior given in Eq. (62). However, one can wonder the next questions: How 
can the free particle know that it is restricted to the interval \x\ > a where the behavior of 
the function is "good"?, How an state of a physical system can be represented by a function 
which is mathematically defined as zero outside some interval? 

At this point, let us recall the meaning of the wave function in Quantum Mechanics. In the 
first place, the wave function represents the physical state of a quantum system. That is to 
say, it represents a combination of the physical properties like energy, momentum, position, 
etcetera, that can be ascribed to the system. In the second place, \ip(x, t)\ 2 dx gives the 
probability that the particle could be found between x and dx. Then, the wave function 
carries the whole information available for the system. For example, a confined particle is 
restricted to have certain eigenf unctions that belong to the domain of the Hamiltonian, and 
certain eigenfunctions that belongs to the domain of the momentum operator [25, 26, 27, 28], 
see also reference [29]. The most severe restriction is / \tp(x)\ 2 dx < in the entire interval. 

Now, we can give a preliminary answer based in a physical insight to the question quoted 
in the previous paragraph: How can the free particle know that it is restricted to the interval 
\x\ > a? By definition, a free particle is a not restricted particle. Therefore, the answer to 
the question is that there is no way that the particle knows that it is restricted to certain 
interval, at least if there is not an infinite well where the particle is confined. That is, only 
in the case of a particle confined in an infinite well we can set the condition ip(x, 0) = 
outside of the well, and the Physics changes from that associated to a free particle to that 
associated to a confined particle [30] . 

As a conclusion of the previous paragraph, we can state that physically the wave function 
i/jh(x, 0) (when it is defined in the entire interval — oo to +oo ) is not valid for the free 
particle. In fact, the answer is related with the differences between Hermitian and self-adjoint 
operators. Mathematically, an operator consists of a prescription of operation together with 
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a Hilbert space subset on where the operator is defined [24, 25, 26, 28]. That is, the functions 
have to belong to the domain of the operator, and if the operator is defined in some interval 
then the set of functions where the operator is defined, i. e. their domain, have to be defined 
in the same interval. Therefore, we think that the above example is mistake because the 
function tp H (x,0) does not belong to the domain of the hamiltonian operator of the free 
particle. 

Let us explain, the difficulties comes from the fact that in Quantum Mechanics the observ- 
able is represented by operators (in a Hilbert space) and the physical states are represented 
by vectors (wave functions). However, the definitions of both operators and vectors given 
in most textbooks of Quantum Mechanics are very weak. The majority of them define an 
operator as an action that changes a vector in another vector, and after that they define 
Hermitian operators as a symmetric operator. However, there is not any mention about the 
domain of the operator and the differences between a self-adjoint and Hermitian operators. 
Because of this weak definition there are many problems or "paradoxes" in the calculations 
of physical properties, see the examples given in references [25, 26, 28]. To handle these 
problems the concept of self-adjoint extension is reviewed in references [25, 26, 28]. Also see 
references [31, 32]. 

Moreover, we adhere the recommendation of Klein [24], Araujo et. al. [25], Bonneau 
et. al. [26] and Gieres [28]: it is necessary to define always the domain of the operators. 
Therefore, in order to droop up all these problems, we think that it is better to define 
operators in analogy with the definition of a function. Additionally, it is necessary to clearly 
state that, as the cited authors have pointed out, unbounded operators cannot be defined 
on all vectors of the Hilbert space. In particular, it has to be stated that in order that a 
function can be valid as an initial state it has to be inside of the domain of the hamiltonian 
H. On the other hand, the only way to set the condition that a wave function is zero outside 
of some interval is to collocate an infinite well in that interval. 

Therefore, the main point in this subsection is that because the function i/jh(x, 0) is not 
square integrable in the interval — oo < x < oo, then it does not belong to the domain of the 
Hamiltonian operator of the free particle. Therefore the function 0) can not represent 

a state of the free particle. This means that the statements 1) and 2) are not true at same 
time. 
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VI. CONCLUSION 



From the work made in this paper, we can conclude that the evolution operator method 
is an efficient method to calculate the evolution of a wave function. This method requires, 
at first instance, the factorization of exponential operators. The factorization allows us to 
apply the exponential operator individually. We have shown how this method works and 
apply it in elementary unidimensional cases. 

As you may guess, all methods have their troubles and limitations. One trouble of 
the evolution operator method is that it is not always possible to find the factorization 
of the exponential operator. Another trouble with this method is that it is not always 
possible to group the evolving function in a single expression as we show in the Appendix 
A. However, this method is increasingly used by many authors. For example, we can recall 
the work of Balasubramanian [12] who discussed the time evolution operator method with 
time dependent Hamiltonians. Also see reference [13]. 
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APPENDIX A 

Here we show another way to factorize the harmonic oscillator. In this case, we define 
the function as: 

F(0 = e^ A+ ^ = e MO<y«)V«f (Al) 
Remember that C = [A, B\. By differentiating Equation (Al) with respect to £, we obtain: 

dF(0 



and 

dF(0 



(A + B)F(0, (A2) 



F(0, (A3) 
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By equating Equations (A2) and (A3), we obtain the following system of differential 
equations: 



dj_ 





= 1, 




= 1, 


d^ 


do dh 
f— H 

J dC di 


= o, 



(A4) 

subjected to the initial condition: 

F(0) = 1. 

By solving Equation (A4), with the initial condition, we obtain the following solutions: 



/(0 " ^sl2^(-fj^y (A5) 

9(0 = -^tanh(^), (A6) 

h (0 = ^ In [sec/* 2 (£V7)] , (A7) 
and we obtain the factorization we were looking for. 



1. Aplication 

The trouble with the factorization given in Eq. (Al) is that at some time we have to 
apply the exponential, which contains the operator C, of the form: 
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with ai a constant. The application of this exponential means to apply(xp) n {e x } which 
produces the following set of polynomials A n (x) (n = 0, 1, 2, . . . ): 

A = 1 

A 1 = 2x 2 

A 2 = Ax A -Ax 2 

A 3 = 8x 6 - 24x 4 + 8x 2 

A 4 = lQx 8 -9Qx 6 - 112a; 4 + lQx 2 (A9) 

However, we were not able to find the generating function. 
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